A Majorana Fermion t — J Model in One Dimension 

Diptiman Sen [| 

Centre for Theoretical Studies, Indian Institute of Science, 
Bangalore 560012, India 

and 

B. Sriram Shastry 
Physics Department, Indian Institute of Science, 
Bangalore 560012, India 

Abstract 

We study a rotation invariant Majorana fermion model in one dimension 
using diagrammatic perturbation theory and numerical diagonalization of 
small systems. The model is inspired by a Majorana representation of the 
antiferromagnetic spin-| chain, and it is similar in form to the t — J model of 
electrons, except that the Majorana fermions carry spin-1 and Z2 charge. We 
discuss the implications of our results for the low-energy excitations of the 
spin-i chain. We also discuss a generalization of our model from 3 species of 
Majorana fermions to species; the SO {A) symmetric model is particularly 
interesting. 

PACS numbers: 71.10.Fd, 75.10.Jm, 75.50.Ee 



^E-mail address: diptiman@cts.iisc.ernet.in 
^E-mail address: bss@physics.iisc.ernet.in 



1 



I. INTRODUCTION 



In a recent paper [|I|, we used a representation of spin-| in terms of three 
species of Majorana fermions 0, ^ in order to study the antiferromagnetic 
spin-i chain. The Majorana represenation has an advantage over other rep- 
resentations (such as the Schwinger boson or fermion representations [Q, |^) 
in that one does not have to impose a constraint on the total particle num- 
ber at each site (see however ref. 0). It is also rotation invariant unlike the 
" drone fermion" and the Holstein-PrimakofF boson representations [0, ^ . 

For the spin-| chain with isotropic nearest neighbor interactions, the 
Majorana representation followed by a rotation invariant Hartree-Fock (H-F) 
analysis [|T| leads to a picture of the low-energy excitations of the spin-i chain 
which is qualitatively similar to that obtained by other methods [P, |l^, |lT|. In 
particular, we find that the excitations are described by a two-parameter con- 
tinuum in the (g, u) space; for each momentum q, the low-energy spectrum 
has a range of energies u as if the excitations are made up of two particles 
(called "spinous"). We also get reasonable dynamic structure functions and 
susceptibilites at all temperatures if we introduce some phenomenological 
structure functions. We should note however that our Majorana fermions 
carry spin-1 unlike the "standard" spinous with spin-|. 

The positive features of the Majorana representation encourages us to 
study the fluctuations about the H-F state of the spin-| chain. More gen- 
erally, it seems to be interesting to examine a strongly correlated Majorana 
fermion model in one dimension and contrast its properties with the much 
better studied electronic systems like the Hubbard model. Such an analysis 
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would also be useful for other possible applications of Majorana fermions such 
as the Kondo problem 0]. In this paper, we therefore study the t — J model 
with Majorana fermions; the electronic version of this model has played a ma- 
jor role in theories of strongly correlated systems like the high-temperature 
superconductors. 

An outline of our paper is as follows. In section II, we briefly recall the 
Majorana representation of spin-| and the H-F analysis of the antiferromag- 
netic chain given in our earlier paper . This motivates a study of the t — J 
model which is introduced in section III. We present the Feynman rules for 
the propagator and the vertex, and compute the one-loop correction to the 
propagator. In section IV. A, we compute the two-loop correction to the prop- 
agator; we find the remarkable result that the on-shell correction is of the 
same form as the tree level dispersion relation. In section IV.B, we compute 
the two-loop correction to the dynamic structure function. The result can be 
used to perturbatively improve the power law of the equal-time correlation 
function and the ground state energy of the spin-| chain from the values 
obtained at the H-F level. In section IV. C, we study the one-loop correction 
to the vertex. In section V, we discuss the symmetries of the t — J model and 
numerically analyze the spectrum of small systems using exact diagonaliza- 
tion. In section VI, we generalize our model from 5*0(3) to SO{N), and we 
briefly examine the 5*0(4) case which is particularly interesting. Finally, in 
section VII, we summarize our understanding of the t — J model. 
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II. MAJORANA FERMIONS AND THE 
ANTIFERROMAGNETIC SPIN-i CHAIN 



At each site n, the spin operators Sn = ct„/2 can be written in terms of 
the Majorana operators 0„ as [|l], ^ Q 

and < = - z 0^ 0^ . (1) 

(We set Planck's constant equal to 1). The hermitian operators 0^ (with 
a = x,y, z) satisfy the anticommutation relations 

{ 0^ , 0^ } = 2 5mn 5ab • (2) 

Note that there is a local Z2 gauge invariance since changing the sign of 
0n does not affect Sn- We will therefore say that 0n (or any odd power 
of it) carries a Z2 charge. Let us define the trilinear and hermitian object 
ipn = -i4>n<Pl4>n- Then [a^,ipn] = 0, and {ipm,i'n} = '^Smn- Under rotations, 
(pn and (Tn transform like vectors (spin-1 objects), while ipn remains invariant. 
On the other hand, ipn carries a Z2 charge while ct„ is Z2 neutral. Thus we 
have two different composite operators, a„ and ipn, which carry spin and 
charge respectively. 

For a system with L sites, it is known that the minimum possible di- 
mension which allows a representation of the form given in equations (|I]-0) 
is 2'^+'^/^!, where [L/2] denotes the largest integer less than or equal to L/2. 
For L sites with a spin-^ object at each site, the Hilbert space clearly has a 
dimension 2^. Thus the Majorana representation of spin-| objects requires 
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us to enlarge the space of states; the complete Hilbert space of states is 
given by a direct product of a 'physical' space and an 'unphysical' one. The 
operators (?„ act only on the physical states, while the 0„ mix up different 
unphysical states. 

We now consider the Heisenberg antiferromagnetic chain with the Hamil- 
tonian 

H = J ^ Sn ■ Sn+i , (3) 

n 

where J > 0. We use periodic boundary conditions S^+i = .Si. The spectrum 
of @j is exactly solvable by the Bethe ansatz; the ground state energy per 
site for large L is given by Eq/L = (— ln2 + 1/4)J = — 0.4431J. The lowest 
excitations are known to be four-fold degenerate consisting of a triplet {S = 
1) and a singlet {S = 0) [|lOl- The excitation spectrum is described by a 
two-parameter continuum in the {q, u) space, where — vr < g < vr. The 
lower boundary of the continuum is described by the des Cloiseaux-Pearson 
relation 

i^iiQ) = ^ I sin g I , (4) 
while the upper boundary is given by 

uju{q) = ttJ I sin I I . (5) 

We can understand this continuum by thinking of these excitations as being 
made up of two spin-| objects ("spinous") with the dispersion [jlO| 

^s{q) = ^ sin g , (6) 

where < g < vr. A triplet (or a singlet) excitation with momentum q is 
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made up of two spinons with momenta qi and q2, such that < qi < q2 < tt, 
q = qi + q2, and uj{q) = uj^iqi) + uJsiq2)- 

The Majorana analysis of this system proceeds as follows We write 
@ in terms of Majorana operators and then perform a H-F decomposition. 
Thus 

H = - J ( 4>l4>Wn+i4>l+i + cycl. perm, (x, y, z) ) 

n 

{<Pn<Pn+i){<Pl<Pl+i) + cycl. perm. {x,y,z) ] . (7) 

For a rotation and translation invariant H-F analysis, we have g = , 
where g has the same value for all n and a = x,y, z. (Our earlier paper [|1| 
follows slightly different conventions) . The Fourier expansion of 0„ is defined 
as 

^ E [ e^'" + ^1. ] , (8) 

0<(?<7r 

where {^ag^^Ig'} = ^ah^qq'- We will work with antiperiodic boundary condi- 
tions for 0^ and even values of L in order to eliminate modes with q equal to 
and TT. In equation d^), g = 27r(p — l/2)/L, with p = 1, 2, L/2. In the 
limit L ^ oo, we get 

^ = E E blqbaq + 3LJ ( ^ - f ) , (9) 

where the Majorana fermions have the dispersion Ug = v sin q, with v = 2gJ. 
The value of g is determined self-consistently tohe g = 2/7t. The H-F ground 
state energy is therefore 

= - ^ J = - 0.3040 J , (10) 
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which is greater than the exact value mentioned above. The "spinon" spec- 
trum has the same form as in (|]), except that we get v = AJ/n instead of 

'Vexact = TtJ/2. 

We can go on to show that the Majorana fermion has spin-1, and a two- 
fermion state therefore has = 0, 1 or 2 in general However the state created 
by Sg = J2n S^e~^'^^ , where < g < vr, has the form 

si\Q) = -I E ^U^L-fc 1 ) ' (11) 

and has S = 1. We thus obtain a two-parameter continuum of triplet exci- 
tations as in equations (j^-]^), with a prefactor 4/7r instead of 7r/2. 
Finally, the equal-time two-spin correlation function is given by 

= ( I 5„ ■ 5o I ) = ^ for n = , 

This does not agree with the correct asymptotic behavior of which is 
known to oscillate as (— l)"/n. In particular, the H-F static structure func- 
tion S{q) = '}2nGn^~"^'^ does not diverge as g — tt in contrast to the cor- 
rect S{q) which has a logarithmic divergence at vr. (Note that we do get 
Z^nG'n = 0, as expected for a singlet ground state). We will show in section 
IV. B that two-loop effects effectively reduce the power governing the asymp- 
totic decay from 2 to 1.75 which is somewhat closer to the correct value of 1. 
At the same time, the ground state energy per site is reduced from — 0.3040J 
to — 0.3338J which is also closer to the Bethe ansatz value of — 0.4431J. 

One can now consider fluctuations about the H-F ground state by doing 
loop calculations. However, instead of studying only the Hamiltonian (|^) as 
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is sufficient for the spin-i chain, it is useful to study a more general model 
which has the same structure but has two parameters instead of one; the 
parameters are a hopping amplitude t and a quartic interaction J. This is 
the subject of the following sections. 

III. THE MAJORANA t - J MODEL 

We consider the Hamiltonian 

—it J 
H = CC+i - T E ( C0n0n+i<^n+i + cycl. pcrm. (x, z) ) , 

a,n n 

(13) 

with t chosen to be positive, and we perform a perturbative expansion with 
the quartic term. To begin the diagrammatic analysis, we generalize the 
Fourier expression to the interaction picture field 
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a ^i(qn-Ljjqt) 



T ^ -rq ^ ' 
-7r<ij<7r 

where (j)^ = bag if < g < vr , 

= bl_g if -7r<g<0, (14) 

with 

ujq = t sing (15) 
for all q. Then we obtain the propagator 

(0| T rSt) 0^,(0) |0) = I C^'iq^t) = I S^" G{q,t) , 



q\ J -r ~q\ 

coo 



/oo 
dt G{q, t) e^"^* 
-oo 



where i] is infinitesimal and positive, and 6{q) = lif0<g<7r and —1 if 
— TT < g < 0. For loop calculations, it is convenient to define a propagator 
even for values of q not lying in the range [— tt, tt]. To do this, we first define a 
momentum q = q + 2n7c where the integer n is chosen such that — tt < q < tt. 
Then we define G{q,Lj) = G{q,Lj) using (|T^). The propagator is shown by a 
solid line in figure 1 (a). 

The vertex shown in figure 1 (b) is obtained by Fourier transforming the 



quartic term in (|l^). The Feynman rule for the vertex is found to be 

iT{ai, gi, cji; aa, q2, Qs, t^s; ^4, Qa, ^^4) 

= ^(27r)2 6p{Y^q,) diY^oo,) 4 J cos( Vg,) ■ 

i i i 

■ ' 5"^"^ 6^^'^' sin ^(gi - q^) sin ^{q^ - q,) 

+ cycl. perm, (aa, ga; as, ^s; «4, ^4) ] , (17) 

where the spin indices ai to 04 can take the values x, y, z, and the momenta 
gi to g4 need not lie in the range [— 7r,7r]. The periodic 5-function is defined 

as 

5p{q) = E ^(9 - 2^^) • (18) 

n=— 00 

The expression in (|l7|) is antisymmetric under the exchange of any two labels 
{ai,qi,u!i) and {aj,qj,ujj); it also vanishes if all the indices Oj are equal. 

We now compute the simplest loop effect, namely, the one- loop contribu- 
tion to the propagator shown in figure 2 (a). It is called one-loop because 
there is one energy-momentum we have to integrate over. To this order in 
J, the self-energy is found to have the energy independent form 

SW(g,cu) = — sing , (19) 
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where the superscript (1) denotes the order of the loop. Thus the dispersion 
relation changes from (|15D to 

4J 

oj = (t H ) sing . (20) 

We will therefore use the expression (|20|) in the propagator ( |T6D for all the 
loop calculations below. Note that we can recover the antiferromagnetic spin- 
i chain by setting t = in (]T3|); equation ( pOD then gives us precisely the 
H-F dispersion discussed in section II. 

IV. LOOP CALCULATIONS 

A. Two-Loop Contribution to Propagator 

We will now compute the two- loop diagram shown in figure 2 (b). The 

two energy integrals can be easily done using the identities 
da; 1 1 i 



oo27iuj — a + irjuj — (3 — ir] (3 — a + i-q 

^ du 1 1 ^ 

-oo 2tt uj — a + i?] ui — (3 + irj ' 

if a and (3 are real. 

We then obtain the following expression for the self-energy 

4J2 r r ,1 ,1 si n^[g + |(/i + /2)] - k)] 

(22) 



where we take the upper sign {irj) in the denominator if — tt < /i, /2 < and 
< /i -|- ^2 + <? < TT, and we take the lower sign {—irj) if < /i, /2 < and 
— TT < /i + /2 + g < 0. It is clear at this point that 

= S(2)(g,^) ; (23) 
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this property of the self-energy can be shown to be true to all orders in J. 

Further, S*^^)(7r — q,uj) = S'^^^(g, c<j). Now let us choose < g < vr and find 

the on-shell dispersion relation to order J^, namely, 

AT 

= + — ) sing + S(2)(g,cu) . (24) 

71 

To this order in J, we can set uj = {t + AJ/n) sing in the second term on the 
right hand side of (|2^) or, equivalently, in the denominator of (^). We then 
find that the denominator in (|2^) never crosses zero in the given ranges of li 
and I2', thus we can drop the ±iri and the integrals are purely real. We then 

numerically find that has the remarkably simple form 

4 J 4 P 

S(2)(g,(t + _)sing) = - 0.467 ——^ sing. (25) 

for all g in the range [0, vr]. Thus the dispersion relation to order is 

4J 

u = (t + 0.189 tt) sing. (26) 

^ TT t + — ^ 

We find it surprising that the form of the dispersion relation remains the 
same even at two-loops, and suspect that this may be true to all orders in J. 

B. Two-Loop Contribution to Dynamic Structure P\inction 

We will compute the two-spin correlation function 

S,,{q,uj) = Fourier transform of (0|5^(t)5^(0)|0) (27) 

to two loops. To any order, we can show that this function remains invariant 
under (g, u) — > (— g, —u). We can obtain the static structure function (equal- 
time correlation function) Szz{(l) by integrating 

/OO ijfj 
Szz{q,u;) e^-\ (28) 
-00 /TT 
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and taking the limit t — > 0+. This is a function of |g|, so it is sufficient to 
compute it for < g < TT. 

The lowest order result for the correlation function is obtained from the 
one-loop diagram in figure 3 (a). After doing the energy integral, we obtain 

1 



00 



I 
2^ 



dli 



1 

cu/j + Uq_i^ -uj + ir] 



if < /i, q-li < TT 

if - TT < Zi, q-h < . 

(29) 



For — 7r < g < TT, we then obtain 



(30) 



The Fourier transform of this gives the spatial correlation function in (p!2D. 

At two loops, we have to compute the diagram given in figure 3 (b). After 
performing the two energy integrations, we arrive at the expression 



where ei 

62 



iJ 
4^ 



TT J —TT 

1 



(ei - iri){e2 - ir]) 
1 

(ei +ir]){e2 + ir]) 
-1 

(ei - iT]){e2 + ir]) 
-1 

(ei +ir]){e2 - ir]) 

OJ + OJli — O^h+q , 



dli dl2 sin ^(/i + h) sin g + ^(/i + h] 



if < li,l2,-{h + q),-{h + q) < TT 



if < -/i, -/2, /i + g , ^2 + g < 
if < /i, -/2, -( /i + g ), /2 + g < TT 

if < -Zi, ^2, ^1 + 9, -(^2 + g) < TT 



2+q • 



(31) 
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We then get, for < g < vr 

J 



2n^it+^) ^^^^ 



9 

Jo 



q COS 
dli dl2 



Wi + k) cos q - + h 



' sin li + sin I2 + sin(g — li) + sin(g — ^2) 
We find analytically that /(g) vanishes as g — and numerically that 

dq I{q) = . 



(32) 



(33) 

These are consistency checks following from the facts that the ground state 
is a singlet and that the two-spin correlation at the same spatial point is 
equal to 3/4; we already know that the one- loop correlation in equation (|l^) 
satisfies these checks. 

We now use equation (^) to derive some interesting numbers relating 
to the antiferromagnetic spin-| chain. First of all, we can show analytically 
that /(g) is finite for all g, while I'{q) diverges logarithmically at g = vr with 
coefficient 1, namely. 



/'(g) = ln|7r — g| + nondivergent terms , 
/(g) = /(tt) + (g — vr) In jvr — g| as g — vr 



(34) 



At long distances, the leading term in the spatial correlation function G. 
3{0\S^S^\0) takes the form 
r'^ dq 



[ — I{q) cos(gn) 

Jo TX 



mn + C(— j as n 



00 



(35) 



After adding this to the one- loop result, we see that the long distance corre- 
lation function has an oscillatory term going as 



Gn 



J 



Inn + 



(36) 
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where the dots indicate contributions from more than two loops. If we now 
assume that these higher order terms come with the right numerical factors 
to turn the sum into an exponential series, we see that the long distance 
correlation decays as (— 1)"?7,~", where the exponent a goes as 

to order J. For the spin-| chain, we must set t = 0; this gives a = 1.75 to 
this order. 

The second interesting number for the spin-| chain which we can derive 
from (^) is the ground state energy per site; this is equal to JGi for t = 0. 
On numerically integrating we find the two- loop result 



= - — / dq I{q) cosq = - 0.0298. (38) 

877"' JO 

On adding this to the one-loop result, we get the value — 0.3338J. 

C. One-Loop Contribution to Vertex 

For completeness, we will mention the one-loop correction to the vertex. 
Let us choose two of the spin indices to be x and two to be y. From (0), the 
zero-loop form of the vertex is given by (27r)^ times the energy-momentum 
conserving (5-functions times 

ir(°) = iU cos{^J2^i) sin^(gi-g2) sin^(g3-g4). (39) 

i 

The one-loop correction iT^^\x, qi,uji; x, q2, (^2; y, gs, co^, y, q^, UJ4) is given by 
the sum of the three diagrams shown in figure 4. On doing the energy integral, 
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we find that the contribution of figure 4 (a) is 

-i8J^ cos(^^gi) sin^(gi-g2) sin^(g3-g4) ■ 

i 

■ r ^ Sin2[/ + l(gi + g2)] ■ 

■ ^ if < /, -(/ + (?1 + ^2) < TT 

■ i if < -/, / + gi + g2 < 71- • (40) 

+ CUl + U;2 - t^i+gi+ga + 

The contribution of figure 4 (b) can be obtained from equation (^OD by chang- 
ing the coefficient 8 to 16 and cyclically replacing q2 q-i (I2 and 
UJ2 —^^^3- The contribution of figure 4 (c) can be obtained from equation 
( ^OD by changing 8 to 16 and replacing q2 —>■ ^3 ^ (I2 and UJ2 ^ uja- 

V. SYMMETRIES AND NUMERICAL RESULTS 

A. Numerical Results 

We can numerically study the spectrum of our model by exact diagonal- 
ization of small systems. To do that, it is useful to know all the symmetries 
of the model. Some of the conserved quantum numbers are the total spin 
5*^ and any one of its components, say, 6*2, the total momentum q modulo 
27r, and parity P = ±1 which arises from the symmetry of the Hamiltonian 
under 

4 ^ (-1)" $L+l-n . (41) 

In addition, there is a Z2 quantum number defined as follows. Consider 
r = ipiip2 ■ ■ ■ if — is even 
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= i ■ ■ - tpL if ^ is odd , (42) 

satisfying T'^ = T^^ = T. This operator anticommutes with each of the (pn 
and therefore commutes with the Hamiltonian (0). Hence the eigenvalue of 
r = ±1 is a good quantum number. We will define F = 1 for the ground 
state of the t — J model; we can ensure this by introducing a — sign in the 
definition ( ^21) if necessary. 



There are a few selection rules and energy degeneracies connecting some 
of these quantum numbers. We will see below that the ground state has 
g = 0, and we choose F = 1. We can now obtain various excited states by 
acting on it with a certain number of Majorana operators as defined in (P). 
Each such operator carries a momentum q which is an odd multiple of vr/L, 
and F = — 1. It is therefore clear that all states must have F = exp(zgL); this 
eiegenvalue is 1 or —1 depending on whether the state has an even or odd 
number of Majorana fermions. Secondly, Majorana operators with momenta 
q and vr — g carry the same energy by parity. Thus states with an odd number 
of Majoranas, i.e. with F = — 1, must have an energy degeneracy between 
total momenta equal to q and it — q. States with F = 1 must have an energy 
degeneracy between momenta q and —q. 



For a numerical study, it is more convenient to rewrite (|T3D in the form 



1 ^ 

if = - ^ ( - it ^pni^n+l + J ) ^n- <?n+l • (43) 
^ n=l 

As mentioned in section 11, we use periodic boundary conditions for an and 
antiperiodic for ip^- We diagonalize if in a basis consisting of a direct product 
of states of the form ® |q;j), such that the operators a„ and ipn act only 
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on and \aj) respectively. In order to study the spectral flow from the 
pure- J model to the pure-t model, we introduce a parameter x lying between 
and 1, such that J = 4(1 — x) and t = 4x. Thus 

N 

H{x) = ^ {I - X - ix ^n^n+l ) " ■ (44) 

n=l 



We have obtained the eigenvalues of (|^) for L = 4 and 6, for 11 equally 
spaced values of x from to 1. All the conserved quantities discussed above 
have discrete eigenvalues; hence these remain invariant as x changes. Nu- 
merically, we only kept track of the eigenvalues of total spin S' = 0, l,...,L/2 
and total = 0; whenever necessary, the eigenvalues q, P and F can be 
deduced by continuity arguments from the exact analytical solutions known 
at X = 1. The energy eigenvalues in each 5* sector are shown in figures 5 
(a-c) for L = 4, and the lowest few eigenvalues in each S sector are shown 
in figures 6 (a-d) for L = 6. We should remark here that the degeneracies of 
the various levels have not been shown, and that we have not distinguished 
between true crossings and avoided crossings in these figures. 

To get a feeling for the elementary excitations, let us discuss six low-lying 
states marked a — / on the figures; these include the three lowest states a, b, c 
with 5 = in figures 5 (a) and 6 (a), the two lowest states d, e with S = 1 
in figures 5 (b) and 6 (b), and the lowest state / with 5 = 2 in figures 5 
(c) and 6 (c). The energy dependence of these six states can be seen to be 
quite similar for L = 4 and 6. The ground state, marked a, is unique for 
all values of x (except x = where it has a degeneracy of 2^/^); it has spin 
S = 0, momentum q = 0, and F = 1. The next two states in the S = 
sector, marked b and c, have F = 1 and —1 with degeneracies of 1 and 2 
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respectively; these two states exhibit a true level crossing between x = and 
1, so that b is lower than c near x — 1 and vice versa near x — 0. The two 
states with S — 1, marked d and e, have F = — 1 and 1 with degeneracies of 
2 and 1 respectively. These also exhibit a true level crossing, with d being 
lower than e near x = 1 and vice versa near x = 0. Finally, the state with 
S = 2 marked / has F = 1 and is nondegenerate. 

The composition of these six states can be easily understood at the non- 
interacting point X — 1. At this point, the ordering of energies is given by 
a<d<b — e — f<c. The ground state a is the empty state. State d 
contains a single Majorana fermion with spin 1, with momentum equal to 
either t:/ L or tt — tt/L; hence the double degeneracy. The state b contains two 
fermions in a spin-0 combination, one with momentum tt/L and the other 
with momentum tt — tt/L; hence the total momentum is tt. States e and / 
have the same composition as b, except that they have spins 1 and 2 respec- 
tively. State c has three fermions in a spin-0 combination, two with momenta 
tt/L and n — n/L, and the third with momentum either Stt/L or tt — 3n/L; 
the double degeneracy is due to the two-fold choice for the third fermion. If 
we now move from x = 1 to x = 0, all these states get "dressed" with an even 
number of fermions. At x = 0, the energy ordering is a = b = c<d — e<f. 

Although the system sizes are not large, we can draw the following qual- 
itative conclusions from these figures. First, the states evolve smoothly from 
X = to X = 1 with no abrupt changes in between. In each spin sector, the 
lowest energy states at a; = are mainly composed of the lowest energy states 
at x = 1, and vice versa. Finally, the complex pattern of level crossings for 
small values of S seems to suggest that the model is nonintegrable for x not 
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equal to or 1. 



B. Conformal Field Theory: A Conjecture 

It would be useful to understand the low-energy excitations of the model 
in terms of conformal field theory; amongst other things, this would lead 
to a simpler derivation of various correlation functions (see ref. and 
references therein). We would like to advance a conjecture in this direction. 
Before doing that, we must consider the two limits of the Hamiltonian (^) 
which are exactly solvable. 

For X = 1, we have three uncoupled Majorana fermions with the same 
dispersion (pTS)). The low-energy excitations (modes with momenta q close 
to or tt) have velocity t = 4 and are governed by a conformal field theory 
which is an SU {2)2 Wess-Zumino-Witten (WZW) model with central charge 
c = 3/2. 

For X = 0, the unphysical states decouple completely. The physical states 
(each of which have an unphysical degeneracy of 2'^/^] due to the spinless 
Majorana field ipn) are solvable by the Bethe ansatz; the low-energy physical 
excitations have the velocity vr J/2 = 27r and are governed by a SU (2)i WZW 
conformal field theory with c = 1. The a; = limit is somewhat singular 
due to the complete decoupling of the unphysical states. Let us therefore 
examine what happens if x is nonzero but small. We can then do degenerate 
perturbation theory to first order in x. For instance, consider perturbation 
theory amongst the 2^^/^' ground states which are degenerate for x = 0; we 
denote these states by the direct product ®\oi), where \l/o is the physical 
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ground state and a can take 2'^/^! values. By rewriting 0„ = Oni^n and using 
the Bethe ansatz value 

e = (^o|^n-c?n+i|^o) = - 1.7726, (45) 

we find that the first term in the Hamiltonian ([T3|) can be written as the 
perturbation 

y = -ixe V'nV'n+l • (46) 

n 

This can be diagonalized by Fourier transforming as 

= M E [4^^'" + C.^"^'"] • (47) 

V ^ 0<g<7T 

Then 

V = - Axe V sin q dcg + . (43) 

0<(?<7r " 

Thus the spinless sector with Z2 charge has low-energy excitations with ve- 
locity —Axe. These are described by a conformal field theory with c = 1/2. 
Thus the spin and charge excitations have completely different velocities. 

The question now is what happens in between the two limits. Although 
our numerical results are limited to L = 4 and 6, they suggest that both the 
spin sector (for instance, states with S > and F = 1) and the charge sector 
(states with 5* = and F = —1) remain gapless for all values of x; there does 
not appear to be a quantum phase transition at any point between x = 
and 1. It is then natural to conjecture that the low-energy sector is generally 
described by the product of two conformal field theories which have different 
velocities; the spin sector by a SU {2)i WZW model with c = 1, and the Z2 
charge sector by a single Majorana fermion with c = 1/2. If this is correct, 
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it would be somewhat reminiscent of the one-dimensional Hubbard model 
away from half-filling; the low-energy excitations of this are governed by the 
product of two conformal field theories which have different velocities if the 
on-site interaction f/ 7^ 0; the spin sector is again described by a SU{2)i 
WZW model while the U{1) charge sector is described by a Gaussian field 



theory with c = 1 [y, [T|] . 

VI. SO{N) t - J MODEL 

It is possible to generalize the t — J model with three species of Majorana 
fermions to a model with species. In terms of an interpolating parameter 
X, we can write a SO{N) symmetric Hamiltonian in the form 

N 

H = - CC+l - (1-^) E E C 0n C+1 0n+l , (49) 

n a=l n l<a<b<N 

where the operators 0^ satisfy the same anticommutation relations as in (|^), 
except that the flavor indices a, b can now take values. The Hilbert space 
for L sites has the dimensionality 2^^/^ if L is even. For x = 1, we have 

noninteracting Majorana fermions with the dispersion Ug = 4 sing; the 
low-energy excitations are therefore described by a c = A^/2 conformal field 
theory. We will now examine two special cases, N = 2 and = 4, for which 
the antiferromagnetic limit x = is also well understood. 

For N = 2, the model is equivalent to the XXZ spin-i chain. This can 
be shown as follows. We first combine two Majorana operators to produce 
an annihilation operator for a spinless Dirac fermion. 

dn = ^ {ct>l + ict>l). (50) 
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These satisfy the anticommutation relation 

{dm , dl} = 5mn ■ (51) 

In terms of these, the Hamiltonian takes the form 

1 1 

H = 2xJ2 dldn+l + )+ 4(1-X) ( didn - ^){ " - ) • 

n n ^ ^ 

(52) 

A Jordan- Wigner transformation from fermions to spin-| operators then pro- 



duces the XXZ Hamiltonian 12 



H = xY.{ «+i + ^>n+i ) + {l-x)Y. «+i • (53) 

n n 

This model is exactly solvable by the Bethe ansatz for all values of x\ it has a 
quantum phase transition at x = 1/2. For 1/2 < x < 1, the model is gapless 
and is described by a c = 1 Gaussian conformal field theory (the symmetry 
is enhanced from f/(l) to SU{2) at x = 1/2). For < x < 1/2, the model is 
gapped and has a Neel ground state with long range order. 

The case = 4 is more interesting. At x = 0, the model is a direct sum 
of two antiferromagnetic spin-1 chains. To show this, let us first define the 
six generators of 5*0(4) at each site. 



K' = \ ■ (54) 

Now we use the homomorphism 5*0(4) ^ 50(3) x 50(3). This can be proved 
by defining the linear combinations 
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(55) 
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These generate two commuting 5*0(3) algebras, namely, 

c 

where = 1,2 label the two algebras, a, 6, c = x,y,z, and e^^"^ = 1. We 
can define total angular momentum operators 

^« = E ^an ; (57) 

n 



these commute with the Hamiltonian p9| ) for all values of x. 

At a single site, the Hilbert space is four-dimensional; the four operators 
(f)"' can be chosen to be the 7 matrices used in Dirac's theory of the electron. 
One can verify that 

3 

Lj = g ( ^ - 01</'2</>3<^4 ) , 

Ll = ^{1 + 01020304 ) . (58) 
o 

It is convenient to choose a representation in which these two operators are 
diagonal in the form of 2 x 2 blocks 

/ 

Li = 

V 

LI - f" ° |. (59) 

\^ 3/4 

Thus the upper two components of the Hilbert space transform as the (|, 0) 
representation of (Li,L2), while the lower two components transform as 
(0, |). We now see that, for x = 0, the Hamiltonian for L sites has the 
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block diagonal form 

( Hi 

H = 



(60) 



v 



H2 

where the Hamiltonians Hi and H2 act on two separate 2^ dimensional 
Hilbert spaces, each corresponding to a spin-| chain. Here 

Ha = L (61) 

n 

for a = 1,2. We already know that this can be solved by the Bethe ansatz; 
the block diagonal form of (|6^) implies that each eigenvalue will have a two- 
fold degeneracy. Thus the 5*0(4) t — J model is exactly solvable at both 
a; = and 1, and one can investigate how the spectrum interpolates between 
the two. We will not pursue this here. 

The Majorana fermions in the SO {A) model carry the spin quantum num- 
bers {Li, L2) = (|, !)• In this respect they may be closer in spirit to the 
Faddeev-Takhtajan spinous (which are spin-| objects) than the Majorana 
fermions in the 5*0(3) model which carry spin-1. To show this more pre- 
cisely, let us define two Dirac fermion operators in the 50(4) model as 



din = ( 4>ln + i(p2n ) , 

(— i)" 

d2n = ^ ( 4>3n + «04n ) • (62) 

We can then verify that the particles created by c/|„ and d^n have the eigen- 
values of the total angular momentum operators {LI, Lg) equal to (|, |) and 
(i, — i) respectively. Thus, for the purely antiferromagnetic model x = 0, a 
fermion operator acting on the ground state of, say, the Li chain will pro- 
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duce states which transform as spin-i under the operators Li; in addition, 
the states will carry a two-fold internal quantum number coming from L2. 

It is interesting to note that the hopping term (proportional to t) in the 
5*0(4) Majorana model is identical to the hopping term in the Hubbard 
model of electrons. However the four-fermion interactions are very different 
in the two models. 

Before ending this section, we would like to mention that a H-F analysis 
of the SO{N) antiferromagnet has been performed in ref. [Q. Their H-F 
decomposition differs from the one we have used in section II. Consequently 
they obtain a much higher value for the ground state energy than us, namely, 
equation ([10|) for = 3, and —JN{N — l)/27r^ in general. 

VII. DISCUSSION 

We have studied a one-dimensional 5*0 (3) invariant t — J model with 
Majorana fermions. At the pure-J end, this describes the nearest neighbor 
antiferromagnetic spin-| chain, while at the pure-t end, we have three non- 
interacting fermions. We have done perturbative calculations to low order in 
the four-fermion interaction. We have also studied the model numerically by 
exact diagonalization of small systems. These studies provide a new perspec- 
tive on the excitations of the spin-| chain by relating it in an "adiabatic" 
and rotationally invariant way to a model of free Majorana fermions. The 
low-energy excitations of the spin-| chain can be thought of as being made 
up of a small number of Majorana fermions. 

The field theoretic description of the low-energy excitations of the model 
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remains unclear. We have suggested that these excitations are governed 
by a product of two conformal field theories which have entirely different 
symmetries. Numerical studies, particularly finite size scaling, of much larger 
systems are required to test this scenario. 

The SO{N) generalization of our model also deserves further study. The 
SO (4) case seems to be specially interesting because it provides yet another 
way of smoothly connecting a model of free fermions to a spin-| antiferro- 
magnet. The 30(4) model may also have applications to the problem of two 
coupled spin-i chains ||T5[ . 
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Figure Captions 

1. The propagator and vertex for the Majorana t — J model. 

2. The one- and two- loop contributions to the propagator. 

3. The one- and two-loop contributions to the two-spin correlation function. 

4. The one-loop contributions to the vertex. 

5. Energies for L = 4. x — Q and 1 denote the pure-J and pure-t models 

respectively. Figures (a-c) show all the energies for total 5* = 0,1,2. The 
curves marked a — f are discussed in the text. 

6. Energy for L — 6. Figures (a-d) show the lowest 20 energies for 5" = 0, 1, 2 
and all the energies for 5" = 3. The curves marked a — f are discussed in the 
text. 
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